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A THEORY OF MEMBRANE PERMEABILITY: I 


INGRAM BLOCH 
THE UNIVERSITY OF CHICAGO 


As a first approximation, a membrane is considered as a potential 
barrier or potential well for the molecules which diffuse through it. 
Making some simple assumptions about the form of that potential bar- 
rier or well, the author sets up and integrates the differential equation 
for the diffusion of the molecules through the membrane under the in- 
fluence of the potential field. 


A plane membrane may be conveniently thought of as a collec- 
tion of molecules, hereinafter known as “membrane molecules”, lo- 
cated between two parallel planes. Associated with this membrane 
and any given substance diffusing through it is a potential function, 
whose value at any point can be taken to be the amount by which the 
potential energy of a molecule of the diffusing substance (the “so- 
lute’) at that point exceeds its potential energy an infinite distance 
away from the membrane. This function may be predominantly either 
negative or positive in the membrane, according as the membrane is 
a potential “well” or “barrier”. 

The normal to the membrane will be defined as the x-axis; the 
flat region occupied by the membrane extends from x = —b tox = +b. 
Symbols to be used are defined below: 


J = grams/cm.? sec. of solute passing through the membrane 
in the + «x direction (J is constant in 2). 


D(a) = diffusion coefficient. 


2 ; CCl —By). 
c(x) = concentration of the solute, in gm./cc. c,=c(Bo). 
de 
oe 
dx 
V(x) = potential function. V(x) =Oifa >B,or“x<—B,. 
dV 
LOE Monet 
dx 


T(x) = temperature. T(B2) EY —T,. 
85 
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k= Boltzmann’s constant = 1.371 X< 10-*° E 
deg.C 


h = membrane permeability. 
4= mean free path of solute molecules. | 
v(x) = average velocity of solute molecules. 
v? (a) = mean square velocity of solute molecules. 
N(x) = number of molecules of solvent (water) per cc. 
NM» («#) = number of membrane molecules per cc. 


Nat) == 0 10k 0, —etuaw nO & 
and b <% <0. 
Nm (x) =N £00) =a 


t,o radius of water molecule (assumed spherical) . 

Tm = radius of membrane molecule (assumed spherical). 
y¥, — radius of solute molecule (assumed spherical). 
Thee Tes 

Te et tant 

m =mass of solute molecule. 


a= volume per cc. occupied by neither water nor membrane 
molecules. 


Equaticn (1) below describes the motion of a diffusing substance 
in a force field. 


OS Mae coi ent ae Nk. (1) 


The solution of this equation can be written 


; x 
C = e-SIV' (a) /kT (2) \de Semi CLV NAEP ay tee (2) 
cima rer aes og 
D(x) 
a ~Bx 


where a is a constant of integration. 


The temperature is proportional to average molecular kinetic 
energy, and may, therefore, be approximately written as 


e Z 
fy = tT; Ai ae V ’ 
(x) a (x) (3) 
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if 


3 
Vj) = 5 kT). 
A more accurate expression for 7 might be 


2 
T(x) =T,-—V(x#t+Az), (4) 
3k 


where A « is inserted because in a liquid the average kinetic energy 
of molecules at x may be influenced by collisions with clumps of mole- 
cules whose centers of gravity are at positions denoted by other val- 
ues of x. The quantity A x is positive or negative according as J is 
negative or positive. 

Such is the approximate nature of this treatment that equation 
(3) will be taken as giving T: 


Vi(ejydx- 1 dV 3) 2 
——————_>- ———= — ~-log (T,-— V (5) 
kT (x) k Pe 2 3k A 
f A ene 
3k 


Equations (2) and (5) yield 


since ¢, = c(—B,) and V(—B,) = 0, 
Ge 


ee (7) 
L212 
Further, 
2 3/2 
a T) — ar V (a) |az 
SS eT I ae ae anne ere a ey (8) 
C1 Co 0 af D (x) 
-By, 
The permeability h is defined by 
h(¢>— G:) =J; (9) 


hence, 
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T 3/2 
ee ee (10) 
2 -3/2 
a= 1b! 
3k 
—__—__—_—__—_—_———— dz 
of D(x) 
An approximate expression for D is 
Doe: (11) 
3 
4 can be approximated as 
1 
,= 7. (12) 
n (Nf? 1 Nm Fra”) 


If the solute is greatly diluted, , may be considered to be influenced 
only by virtue of the membrane molecules’ occupying space which 
would otherwise be occupied by water molecules, or the sum of the 
volumes/cc. occupied by the two types of molecules may be set con- 
stant: 


4 
ate Mm + Ty? Ny) =1 — a; (13) 
this gives 
2 (1 75 a.) Car Nm 
Np = —— sa —; (14) 
4 Nl Ta fig 


Equation (14) will subsequently be used for substitution in equation. 
(12). 


It is approximately true that 
v= a (15) 
strict accuracy would require the use of a multiplicative constant of 
value (determined by the type of molecular velocity distribution) 


about equal to unity. If EH, is the average translational kinetic energy 
of a solute molecule at « = B, or —B,, 


1 lees 
5m pis Bens o- V=skTo—V, (16): 


eS 3kT, — 2V 
v= | ————. (17): 


or 
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Substitution in equation (11) from equations (12), (14), and 
(17) yields 


V 3kT, — 2V 
Dp , (18) 
3m ve Mine ( Roe Ry? Tm3 a 3Ry2 1—«a | 
(mae Ase ee 
whence, 
2 -3/2 
(7 or V ) : = 
3k (3k)3/? 82m 3 = 
—_ = [{ - kT, — V 
D 4 2, 
(19) 
Re? qe 3R,2 1 ere OL 
Rim (Re = ua) SF | 
re An Tape 
Equation (19) can be abbreviated as 
2 -3/2 
Ls el ae ae V 
3k | 
= 2 Saar pae roe ere CNA CSE (20) 
where 
82 Vm (3k)? 
Se (21) 
4 
3S 
a, =- kT, ’ (22) 
2 
hore te 
As = R,? = , (23) 
ie : 
and 
3k ile 
ae oe (24) 
Ane £72 
Equations (10) and (20) may now be combined: 
T-3/2 


can dy [32 (Ge — V(x) )* (Asm (ae) + ay) dx’ 
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or, since %, = N or 0 according as « is between —b and +0 or not, 


= T-*2/ay {04 [3? [2 — V (ae) ]* dc + f® [oe — V (ae)? de] + 
26 
(aN + a4) f?,Ed2 — V(a)]* dey. coe 


It can be seen that the quantities B, and B, determine the 
“spread” of V(x), or the gradualness with which V(x) goes to zero 
outside the membrane. If the force field is symmetric, B, = B2; if 
V(x) = 0 for all |x| appreciable greater than 0, then B, = B, = b, 
and the first two integrals in equation (26) vanish. If the potential 
energy of a solute molecule at distance r from a membrane molecule 


A’ 
is given by A/r®, then V (x) will be roughly rer for% > b +05 
Cas os 
A’ 
d ————— for x < — b — 6, where 6 is a small positive number, 
(2 + b)* 


and A and A’ are constants of the same sign (either positive or nega- 
tive). The larger the quantity s, the smaller the distance B, = Bz 
beyond which V(x) = 0. 

The reader may be dissatisfied with the definition given for h, 
since its magnitude is decreased on account of the first two integrals 
in equation (26), which represent the resistance offered by water 
which is really not part of the membrane at all. However, the values 
of these integrals depend on a property of the membrane; namely, 
the potential function; furthermore, any measurement of membrane 
permeability will perforce include in the permeability the diffusion 
coefficient of a certain amount of liquid on either side of the mem- 
brane.? For application of equation (26) to experimental results, it 
is recommended that B, and B, be taken as the distances from the 
membrane at which c¢, and c,, respectively, are measured. 

Methods of determining V(x) will be discussed in a later paper. 

In conclusion, it seems advisable to evaluate h for a membrane 
whose potential function is constant inside the membrane, zero out- 
side: ” 


an 


V (2) =V, for —b =2e=b,; and 


V(x) =0 for — 0 <a < —band forb<2< o. (27) 


1It is of interest to note from the definitions of a,,4,, and a,, that if B, 
and B, are chosen beyond the points where V(x) > 0 , the first two integrals in 
equation (26) will, where V = 0, be quite independent of the properties of the 
membrane, and may be seen to give the effect of resistance of the solvent alone. 
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In this case, B, = B, = b, and 
T 3/2 


SS, (28) 
a,(a;N + a4) »-2b(a, — V,)~ 
or 
3 2 
28h) =-(; kT “to vs) A 
h=——oO  —_____________, (29) 


aed 1 Tn? dh,” (1 =. a) 
37 b Vm i gs ae N =F 
ee An te. 
As might be expected, h is an increasing function of a, and a de- 
creasing function of m, b, and V,. Equation (29) must be inter- 
preted in light of the fact that 


3 
mpi > Vo; (30) 


this is the condition for validity of equation (3), and must hold if 
equation (3) is not to indicate that negative temperatures exist with- 
in the membrane. According to the present simplified theory, a po- 


3 
tential barrier of height V, > 5 kT, would be completely impermeable; 


potential wells are, ceteris paribus, more permeable than potential 
barriers. 
From equation (29), 


h 3 3/3 = 
a 3kT 3%? € kT, — a) oS (5 kT, — Vo Len. (31) 


OT, 


This expression has the same sign as 
3 ; 
2kT _ 2 kT, Se VY = 4 KT, aP Ve 9 (32) 


which is positive uniess 
—V, 2 43kT,; (33) 


that is, permeability is an increasing function of temperature unless 
the membrane is a potential well of “depth” as great as that indicated 
by expression (33). If inequality (33) is satisfied for some tempera- 
tures, the membrane has minimum permeability when 
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—V,=4kT,, (34) 
since 
h 2 
Sion ie (35) 
ay A. k 
and 
oh 2Vs 
<0 forT,< — ; (36) 
T; k 


Dependence of h on N , the density of membrane molecules, may 
be seen from the following considerations: let ry = 6 Tm. Now, 


Ry? + % ne (7s + Bim)? Tn? 
ke SC ——————— (37) 
Tw? BP 
this quantity is positive if 
B (7s = Ser - (T; ot B as > 0 ’ (38) 
or if 
(p°/? hs Lr; Soe (B°/? pies BP) Tm ° (39) 


Inequality (39) is valid if and only if 6 > 1, or hk is a decreasing 
function of N if and only if Tn < Ty. 


The permeability may be seen to be always positive, since 


AnT,? 


N<1— a, (40) 


by the definition of a. 
The special case of a gas diffusing in the absence of a solvent 


may be treated: if there are no water molecules, expression (40) be- 
comes an equality, and 


3 2 
2 (3k) ~/2 (; kT, — Vo ) Sas’ (41) 


k= — 
31 Vm b R.2 N 

Equation (29) indicates that h is a decreasing function of mem- 
brane thickness and of the mass of a diffusing molecule. 

The quantity 1 — a is the same as the Van der Waalsian b; it is 
a measure of the closeness with which the water and membrane mole- 
cules are packed together. Therefore, h is seen to decrease with in- 
creasing closeness of packing. 

The theory herein presented is, though both inaccurate and in- 
complete, perhaps applicable to some cases of diffusion in biological | 


systems. The author hopes to remove some of its deficiences in later 
papers. 
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OUTLINE OF A MATHEMATICAL THEORY OF THE REMOVAL 
OF MALARIAL PARASITES FROM THE BLOOD STREAM 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


_ By making some plausible assumptions, a set of differential equa- 
tions is proposed which describes the kinetics of interaction between peri- 
odically multiplying parasites and continuously produced phagocytes. 
These equations establish some general conditions which are necessary 
in order that the number of parasites should increase or decrease. Some 
consequences of the equations are discussed, and some experimental pro- 
cedures for the determination of the parameters entering those equa- 
tions are suggested. 


There are certain types of diseases, of which perhaps malaria 
is a typical representative, which are due to infection by parasites, 
those parasites multiplying at rather regular intervals. This regu- 
larity of their multiplication impresses a “periodic” character on the 
whole course of the disease. The purpose of this paper is to outline 
an approach to a mathematical theory of the kinetics of interaction 
between the parasites and the host in this type of disease. While the 
following outline in its present form takes care only of the simplest 
aspects of the phenomena, yet, as will be seen, even in this simple 
form it suggests some new approaches in the experimental study. 
Further elaboration of the theory is likely to bring clearer insight 
into the observed phenomena, and to indicate new ways of experi- 
mentation in this field. 

In malaria, each asexual parasite which lives in a red blood cell 
divides at almost regular intervals ;, called the cycle time, into a num- 
ber M of new parasites. M is known as the merozoite number. Upon 
segmentation of the mature parasite, each new parasite enters a 
new erythrocyte, at least at not too late stages of the disease. Con- 
fining ourselves for the present to stages which are not too advanced, 
we may put the total number N of parasites equal to the number of 
infected erythrocytes. Moreover, the all important feature of the 
phenomena observed is that all parasites of a brood divide almost 
synchronously. Therefore, as a rather good first approximation, we 
may consider that the total number N of parasites in the organism 
increases suddenly at regular intervals of duration +. During each 
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of such discontinuous increases, the total number N suddenly increases 
M times. 

Actually, however, the synchronism of division of the parasites 
is not perfect. There is a certain scattering around the end points of 
the intervals +. Inasmuch as the rate of multiplication p(t) per 
parasite is for all practical purposes a periodic function of period 7, 
it may be represented by a Fourier series: 


2 224% 2ntt 
r= (a sin + 6; cos ; (1) 
rt) T cis 

The function p(t) is practically zero everywhere, except near the end 
points of the intervals 7. We shall subsequently discuss how the co- 
efficients a; and b; may be determined from observations. 

In the absence of any other phenomena affecting the number N , 
such as phagocytosis or any other form of destruction of parasites, 
p(t) represents by definition the quantity . 


(1/N) (dN/dt) =d(log N)/dt. 


Since after a lapse of time 7 the number N increases M times, there- 
fore for any value of t we have 


{ Pp yOr ona (2) 


Due to phagocytosis, the number N of parasites in the intervals 
between the periods of almost synchronous division does not remain 
constant, but decreases. Since phagocytosis of parasitized red cells 
occurs only in special tissues, which in the mammal are located in 
the spleen, the liver, and the bone marrow, the rate of removal of the 
total number N of parasites in the blood may be assumed to be propor- 
tional to the fraction f of the total circulation time which the blood 
takes to go through the phagocytic tissues. If ¢, denotes the time it 
takes the blood to flow through the phagocytic tissues, and T is the 
total circulation time, then 


f=t,/T. (3) 


This simple definition can, however, be used only if we consider 
the whole circulatory system as a simple single circuit, with all parts 
connected in series. Actually, the circulatory system consists of a 
number of parallel circuits, each having a different circulation time 
I; . This requires certain modification of the definition (3) as well as 
of the following deductions. In a subsequent paper it will be shown 
that those modifications do not alter fundamentally the main results. 
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At present we shall for simplicity consider the circulatory system as 
a simple circuit. 

The rate of removal of the total number N of parasites will fur- 
thermore be proportional to the total number n of phagocytes pres- 
ent in all the phagocytic tissues. It will also be proportional to N it- 
self. Denoting by o a coefficient of proportionality, we may set the 
rate of decrease of N as equal to —afnN. The coefficient « is essen- 
tially the probability of the removal of a particular parasitized ery- 
throcyte during a unit of time by a particular phagocyte, when the 
erythrocyte is in the phagocytic tissue. 

We-may now write for the total rate of change of N: 


one t)N N 
pie A) afnN. (4) 


The three quantities a, f , and n are, however, not absolute con- 
stants, but in general depend themselves on the number N. As a re- 
sult of their metabolic activity, the parasites produce various sub- 
stances which in general may affect the above three quantities. We 
shall first discuss the dependence of n upon N. 

In a healthy non-infected organism the number 7 of phagocytes 
is likely to be determined by the resultant of two reactions; a produc- 
tion of phagocytes by the tissue at a rate Q, and their natural “death 
rate” or destruction rate which is proportional to their number. As- 
suming for simplicity that Q is constant, we have, denoting by b a 
coefficient of proportionality, 

dn 


—=Q-— bn. (5) 
dt 


The number n tends asymptotically to an equilibrium value 
Mo =Q/d. (6) 


This value n, may be assumed to be maintained in a non-infected 
organism. 

The metabolites produced by the parasites may, however, in gen- 
eral affect both the rate of production and the rate of destruction of 
the phagocytes. Although very little is known of the latter effect, it 
is definitely known that the presence of parasites stimulates the pro- 
duction of phagocytes. (Taliaferro and Mulligan, 1937). We shall 
therefore consider in this paper only this stimulating effect. In a fur- 
ther refinement of the theory, it may be interesting to consider an 
increase in the rate of destruction of ~ proportional to the rate of 


destruction of NV. 
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The stimulating effect of N upon the rate of production of n will 
be represented by a function 6(N). In general, the effect of N upon 
dn/dt at the time t will be determined not by the instantaneous value 
of N, but by the value which N had at some previous time ¢ — ts 
For the metabolites produced by the parasites may affect not the last 
stages of formation of a phagocyte, but some earlier stages. 

Denoting by A a constant, we may now write 


= 9+ AdIN(t— 4) bn, | (7) 


If a and f were constant, then the system of equations (4) and 
(7) would completely determine N and n at any time, if initial con- 
ditions were prescribed. 

The quantity f however, is, itself likely to be affected by N . The 
metabolites produced by the parasites affect first of all the general 
circulation time due to their effect on the cardiac output. Such an 
effect would change ¢, and T in equation (38) in the same ratio, and 
thus leave f unchanged. However, the metabolites may in general 
produce different vaso-motor effects in different tissues, changing the 
effective cross-section of the capillaries and thus affecting Ut, differ- 
ently than T. From considerations of homeostasis, it is plausible to 
assume that f increases with N. Thus we have 


f=f(N); df/dN>0. (8) 


Finally, we consider the variation of a with N. The quantity a 
depends primarily upon surface factors which determine the adhe- 
siveness of the parasitized erythrocytes to the phagocytes, and upon 
the saturation of the phagocytes, i.e., upon the number of parasites 
contained in each phagocyte. The greater the saturation, the smaller 
a. If S is the saturation, then we may assume that a(S) decreases 
monotonically from a(0) to zero for some large value S,, of S, re- 
maining zero for all larger values. A fair approximation which is 
analytically easy to handle may be obtained by assuming that a be- 
comes zero only for infinite S , but drops rapidly with S . The surface 
factors will be affected by the metabolites, and therefore will be func- 
tions of N. We may put approximately, with § as a constant, 


a(N,S) = a (N) e88, (9) 


The saturation § is, however, itself a function of N , as is seen 
from the following consideration. 
; The total number of phagocytized parasites is equal to nS. It 
mcreases at the same rate as N decreases, that is, at a rate anfN . 


H. D. LANDAHL 97 


But it decreases due to two factors; first, due to the destruction of 
the phagocytes which takes place at a rate bn, and second, due to 
the digestion of parasites by the phagocytes which takes place at a 
rate ynS, n being a coefficient of pr oportionality. Putting 


b+n=C, (10) 
we have 
d(nS) ae dn + dS 
ptrae = a Ns Ca)e (11) 


If f(N) and a (N) as wel! as p(t) are known, then equations 
(4), (7), and (11) determine the quantities N, n, and S as func- 
tions of time. If the solutions of those equations show that N and S 
decrease sufficiently rapidly to zero while n decreases to 1, , this rep- 
resents the phenomenon of recovery from the disease. If on the con- 
trary N increases indefinitely, this represents an eventually lethal 
outcome. Thus the functions f(N), «(N), 6(N), and p(t) together 
with the initial value of N at the moment of infection and with the 
constants A, 6, 6, and Q, the last two being characteristic of the 
organism, determine in principle the course and outcome of the dis- 
ease. 

In the future development of the theory the functions p(t), a(), 
f(m), and.6(N) should be derived from biophysical and biochemical 
considerations. Then, solving the system of equations (4), (7), and 
(11) by suitable approximation methods, we should be able to pre- 
dict the course of the disease under given conditions, and check the 
theoretical derivation of the above four functions indirectly by the 
comparison of the predicted course with the observed. This, how- 
ever, is yet a rather remote possibility. We shall, therefore, discuss 
the much more modest problem of how some of the above functions 
can be determined experimentally by the use of the general theory 
outlined above. Then we shall briefly discuss a general conclusion 
from the above theory in its relation to actual facts. 

Determination of p(t). If the division of the parasites were per- 
fectly synchronous, then the variation of N as a function of time ¢ 
would be represented by the curve A'ABB’ of Figure 1. Actually, due 
to imperfect synchronization, it is represented by a curve somewhat 
like A’A”OB’B. The curves A’A and BB’ represent the solution of 
equation (4) for p(t) = 0. Since far enough from the end of the’ 
interval +, p(t) is practically zero, the actual curve coincides with 
A’A and BB’ in the region AA” and B"B. The curves A”A and BB" 
can be obtained by extrapolation of A’A” and B’B’ until they intersect 
the vertical line AB, which is adjusted so that AO == OB. Let z(t) 
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FIGURE 1 


(Figure 1) denote the absolute value of the vertical distance between 
A”A and A”O, or BB" and OB”. Then p(t) = dz/dt, for p(t) meas- 
ures the deviation of the actual dN/dt from the one that would obtain 
if p(t) = 0. Thus we can graphically construct p(t), and then de- 
termine by standard methods the coefficients a; and b; of the Fourier 
expansion (1). 

Determination of a(N) and n. A direct determination of 
apparently presents great difficulties, and very little is known quan- 
titatively about its variation with time. However, there is evidence 
that the rate of increase of m with time is much slower than the rate 
of increase of N (Knisely, 1944). It is plausible to assume that- 
ao(N) and f(N) are also relatively slowly varying functions of the 
form, 


on (N) = gy oN) da; / aN 0% 
f=fotfi(N); df,/dN > 0. (12) 


Under these conditions we may make a rough estimate of the 
average value afn of afn for the period of one asexual cycle. 

Let N, represent the peak value of N at the end of an interval 7 ,. 
this peak value being graphically corrected for perfect synchroniza- 
tion by extrapolating the line B’B" in Figure 1. In other words, let. 
N be the ordinate of the point B. Let N. be a similar quantity for 
the next interval. During the interval +, the number N, has, accord- 


ing to equation (4), dropped to N,e-*/" and then suddenly increased. 
M times. Hence 


N.= MN,e-"", (13): 
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l njee 
ee) 
afn = ——__—_—__, (14) 


T 


Or 


Thus the knowledge of M, 7, and of two consecutive corrected peak 
values of N gives us afn. For a very rough estimate, actual uncor- 
rected peak values may be used. 

If a, f, and n were constant; then N./N, would be constant for 
any two consecutive peaks. By observing the variation of N./N, as 
a function of the varying average N in the course of the disease, we 
may make a rough estimate of the variation of afn , and hence of afn ; 
with N. 

Determination of a(N). The quantity a is the probability that a 
parasite will be removed by a particular phagocyte in a unit of time. 
Hence ant, is the probability that a parasite will be removed by an 
organ containing n phagocytes. Suppose we inject in a non-infected 
organism a number N, of parasites into the artery supplying a phago- 
cytic tissue, say the splenic artery, and find that N, parasites are re- 
covered in the blood of the corresponding vein. Then 


IN, = IMs 
SNOUT ete cere (15) 
or 
jee (16) 
Qa = 
nt,N« 


If n can be determined by direct observation, then, t, being known, 
we can determine the values of a for different values of N. If, how- 
ever, as remarked above, » varies rather slowly with time, then in 
experiments of that kind n will be practically constant, having the 
value n,. Even if the latter is not known, in this way we can deter- 
mine a(N) up to a constant factor of proportionality. 

A different procedure may also be used, especially for small val- 
ues of N. Suppose that it has been determined that in a certain frac- 
tion of cases, say 50%, the disease develops when a certain number 
N, of parasites is injected into the vein of a phagocytic organ of a 
healthy animal. Suppose now that various numbers N, of parasites 
be injected into the corresponding artery, and it is found that a cer- 
tain value N* > N, will produce the disease in 50% of all cases. This 
means that of N* parasites, N, escaped into the vein, while N* — N, 
were destroyed in the tissue. Hence we have 
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N* = Nj 
nt,N* : 


This gives the value of a(N) for N = N*. If we use different values. 
of N; and N*, corresponding to different percentages of incidence of 
the disease after injection, we may thus gain an idea about the form 
of the function a(N). 
The Determination of n(t). The procedure of injecting N. para- 
sites into an artery of a phagocytic organ and of determining the 
number N,, left in the corresponding vein may be used not on a healthy 
animal, but on an animal with an advanced stage of the disease when 
a certain number N of parasites is already present in the blood 
stream. Such determination may be made at definite times ¢, from 
the initial infection of the animal. The total number N(¢) of the 
parasites already present in the blood stream may also be determined. 
Let the same amount N, always be injected, and denote by a prime 
the amount recovered in the vein of an already infected animal. Simi- 
larly to expression (16) we now have 
N, ee N Bs 
a (N)n(t) = ———. (18) 


a “Pp 


(17) 


where N is the number of parasites just previous to the injection. 
Or, dividing (18) by (16), we have 


GCN) nye Nee IN 
a(Na) No iy aA Noo 


If from experiments on non-infected animals we know a as a func- 
tion of N, then, since in the experiments with already infected ani- 
mals, N is known for any t, relations (18) and (19) enable us to de- 
termine n(t). If direct observations of n are possible, the theory can 
be checked by comparing the directly observed values of n(t) with 
those found by means of equations (16), (18), and (19). 

If the function N(t) is determined by direct count of parasites, 
and thus is known, then equation (7) may be integrated, giving, 


(19) 


‘tty 
nan +e | et AO(N)adt. , (20) 
With known n(t) and N(t), either equation (7) or equation (20) 
may then be used for the determination of the function 6(N 2 
While all the above suggested procedures may meet with certain 
practical and technical difficulties, yet it is important that in principle 
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even this outline of a theory suggests experimental methods for the 
determination of the important functions a(N), f(N), n(t), and 
6(N), which largely determine the course of the disease. Once they 
are determined, we may study effects upon those functions of differ- 
ent physical and chemical agents. 

Some General Consequences. The periodicity of p(t) gives an 
“oscillatory” character to the whole course of the disease, the value 
of N increasing sharply at regular intervals, and then dropping again 
more slowly until the next increase. The final outcome is largely 
(though not solely) dependent on whether the general over-all trend 
of N is to increase or to decrease. To study this general trend, we 
may replace p(t) by an average constant value p , obtained by averag- 
ing p(t) over the interval 7. We shall again consider the case that 
f and nm are varying relatively slowly, so that over not too long a 
period of time they may be replaced by their averages f and mn. Then, 
putting 


y=logN, (21) 
we obtain instead of equation (4), 


Oa im seed 
—=p-—aly)fn. (22) 
dt 


Considering the case that the saturation S is zero, we may put as a 
first approximation, 


a(Y) =O + ay’ Y. (23) 
Integrating equation (22), we find, denoting by y, the value of y for 
t= 0: 
P—a'ofa powin\ 
———— + (v. a ——— ee, (24) 
ina’, fnay 


If into an organism initially free of parasites, in which therefore 
S = 0, we inject at the time t = 0, No parasites, then, if Ny < Ne = 
(p — a'fn) /fno’, and p > a'fn, the number of parasites will asymp- 
totically increase towards the equilibrium value N.. If No > N-, the 
number will decrease towards N,. Due to the variability of f and 7, 
the equilibrium value N, is not constant, but itself varies slowly with 
time. (Figure 2, full lines.) 

If, however, p < a'ofn, then y (or N) will never increase. (Fig- 
ure 2, broken lines.) Hence, a necessary condition for a disease to 
occur and continue, in terms of the present simplified theory is that 
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FIGURE 2 


» > a'ofn. Since, as we have said above, at the beginning of the dis- 
ease we have with good approximation, f = fo; n = m, hence the 
above condition becomes » > a'ofom. This suggests that if p could 
be reduced sufficiently, or if a and 7 could be increased sufficiently 
for a given individual, then we would say that the individual is im- 
mune. Note that p can be reduced by reducing the merozoite number 
M or by increasing 7. 

Now let us consider the effects of the saturation S. If n varies 
very slowly, we can neglect in equation (11) the term with dn/dt. 
Introducing also the expression (9) with (11), we now find 


d & 
= = ay (N) fNe#§ — CS. (25) 
U 


It is readily seen that S approaches in equilibrium value S., 
which is the root of the equation 


a, (N) fNe-8s —CS=0. (26) 


This equilibrium is stable because for S < S, the expression (26) is 
positive, while for S > S, it is negative. Moreover, since a,(N) in- 
creases with N, it can be readily seen graphically that S, increases 
with N also. 

Consider now a hypothetical case in which the number of para- 
sites in an organism has somehow been maintained at an approxi- 
mately constant level, N , so that S has reached a constant level, and 
assume that the parameters in equation (25) are such that S also 
varies more slowly than N when the latter is suddenly changed. If 
We now inject a given number N’ of parasites, it will vary according 
to equation (22) into which, however, equation (9) must now be in- 
troduced. Considering that expression (23) now stands for a(N), 
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we now have: 


dy Se 
ae = (a'o + a'sy) fne*s , (27) 


or, integrated 


Oo. Oo” ne 
= — ae (. + = erm BS ‘ (28) 


Oy 1 


Equation (28) shows that if we start with N =0, S=0 , we 
shall find that after injection the number of parasites will decrease 
relatively rapidly, (Figure 3, solid line) since e-6* = 1 and the coeffici- 
ent of ¢ in the exponent is large. If we start with larger values of N so 
that S > 0, then y (or N) will decrease after the injection much more 
slowly. (Figure 3, dash line.) For still larger N , we shall have a still 
slower decrease. (Figure 3, broken line.) 


FIGURE 3 


From such experiments it should be possible to find something 
about the rate of phagocytosis as a function of the saturation S. We 
should expect a high rate of removal after injection in a healthy or- 
ganism, which should not be the case with a sick organism. 

The above considerations suggest an explanation for the apparent 
high rate of removal of parasites in a newly infected host and the 
low rate in a diseased host necessary to maintain the disease. 

The author is indebted to Professor M. Knisely for a discussion 

aper and for valuable suggestions. 
rae oe was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abott Memorial Fund of the University of Chicago. 
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According to the theory of Boyle and Conway, chloride leaves the 
extracellular fluid to enter cells, particularly muscle, when the plasma 
potassium concentration K rises. Simplified equations are presented to 
describe in terms of K the distribution of chloride expected when the 
quantity of the salt Kp in the body fluids changes in health and disease. 


Notation. 

Volumes are given in liters of water according to the phase oc- 
cupied: J, intracellular; H, extracellular; W=IJ + E. 

Total content of ions in both phases of animal or tissue, in milli- 
equivalents, is designated by boldfaced type, as Cl. 

Concentrations of diffusible ions in milliequivalents per liter of 
water in the phases are expressed in t¢alics, intracellular values being 
marked by the subscript 1, as Cl,. It is assumed that concentrations 
represent activities. 

Symbols for concentrations and contents. 

K = potassium: extracellular, K; intracellular, K,("K,, if asso- 
ciated with Y-; ‘K, , if associated with diffusible anions) . 

N = those non-diffusible cations in H associated with diffusible 


anions. 

Cl = chloride. 

= diffusible anions of any species univalent, in quantity set ar- 
bitrarily equal to K W, or to'Kil+KE. 

a=sum of diffusible anions other than C/ and p. 

Gl eG 

ad=Cl+a+>p. 

Constants, since E and I, Cl and a are held constant. 

Forn: A=n/E; C=n/I; L=n/W. 
For Cl: B=CI/E; D=Cl/I; M=Cl/W. 

Theory. Many features of the Boyle and Conway theory of the 
distribution of electrolytes in muscle are illustrated in Figure 1. The 
organic radical Y- and the inorganic N* are by some mechanism re- 
stricted to the phase I or E indicated. N is mostly sodium; part of Y 
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may include phosphoric acid esters. These ions attract diffusible ions 
of opposite charge: N attracts Cl and other anions a and p; Y at- 
tracts mostly K , most of which may be ionized. These complexes ac- 
count for the effective osmotic pressures of the respective phases. 
These pressures remain equal since the surface S exerts no tension. 


S 
(J) (E) 
Cl.- Gis 
Ve ee (1 Pomel « ep K* a SN? 
Pr p- 
FIGURE 1 


A unique feature of the theory is that the anions Cl, a, and p of 
the cell water exist at a concentration such that their sum approxi- 
mates the value of K in the extracellular fluid. It is assumed that the 
potassium salt of these anions spreads throughout the free water of 
I and E with concentrations in the two phases equal, 'K, being equal 
to K as shown in Figure 1. 

These potassium salts have no osmotic control over body water. 
They diffuse passively along with water as it is moved by other effec- 
tive osmotic forces: into the cells in water intoxication, out of the 
cells in the dehydration of thirst. In dehydration much of the potas- 
sium lost to the urine in excess of that accompanying the nitrogen 
derived from tissue destruction is probably from this source (see data 
of Elkinton and Winkler, 1944). 

Objection has been raised to studying these ionic equilibria as 
did Boyle and Conway in muscle immersed in artificial and unphysi- 
ological media. Yet whereas in such media the equations can be sim- 
plified because the external concentrations can be kept constant 
through the buffering capacity of an infinite volume, application to 
live animals would seem complicated by interacting changes in con- 
centration between phases E' and 1. So far from this being a handi- 
cap, it is only by meeting it, as in the equations to be presented, that 
we realize the possibility of testing the ideas of Boyle and Conway 
through changes in plasma concentrations and of thus avoiding the 
difficulty of measuring cell concentrations through such indirect meth- 
ods as the inulin-space concept. 

The Relation of Cl, to K. The fundamental equation underlying 
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these developments can be derived in at least two ways. The one, a 
special case of the Donnan equilibrium for a membrane with ion 
pores, follows directly from Boyle and Conway. From their general 
equations (5) and (6), K,Cl, = K Cl and Kid, = K d. Division of 
one equation by the other yields Cl,/d, = Cl/d . And substitution from 
their special equation (23), K = d, , leaves 


Ch= (Ci/d) K; (1) 


Thus for equation (1) to obey Donnan concepts not only must 
K = d,, but also must K, =d=N + K. This is true in immersed 
frog muscle over a wide range of changes of concentrations of K in 
the external medium (Boyle and Conway, Table 7). The analytical 
values for A, in fresh muscle from laboratory animals as published 
over the last decade approximate this value. This implies that the 
average valence of Y is one. It is calculated to be 1.09 for the frog 
and 0.805 for the dog (Dean, 1941; pp. 336 and 344). Equation (17) 
from Dean, d, = K(m + 1)/2 m, m being the valence of Y, is the 
exact relation when the mol fraction of K in the external medium is 
small as it must be because of the limitations of toxicity. Equation 
(1) fits the white rat with a high degree of statistical significance 
(Wilde, 1943) even though K, exceeds d markedly. 

A second type of derivation for equation (1) is a natural exten- 
sion of the “‘visual’’ approach presented under Theory. The distribu- 
tion of the diffusible anions in the two phases is predicted by statis- 
tical mechanics. The probability Cl,/d, that a chloride ion will be 
associated with a given 'K, ion in the cell is equal to the probability 
of the analogous association outside the cell between chloride and 
(N + K). Since the concentration of all diffusible anions inside is 
equivalent to 'K, = K, we write Cl, = (Cl/d)K. If the activity co- 
efficients of the anions differ considerably inside and outside the cell, 
the use of a ratio Cl/d in which chloride is predominant in numerator 
and denominator minimizes the error that such differences of activity 
impose upon prediction. | 

An equation so derived is more general in that it applies to situa- 
tions in addition to the special Donnan case described. It may, for in- 
stance, fit some of the theories reviewed by Steinbach (1940, p. 251). 
However, since these have their difficulties, particularly in explaining 
osmotic relations, further discussion seems futile in a paper of this 
iength. These developments are presented as a guide for testing the 
relation between plasma potassium and intracellular or extracellular 


chloride. 
Algebraic manipulation proves the identity, Cl/d = (Cl E + 
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Cl,1)/(d E + dt) =Cl/d. Substitute Cl/d, using d= n + p, into 
equation (1), which when multiplied by W/W gives 


ae Cl/W 
" n/W + p/W 


Consider the special system in which during the development and 
growth of the animal, or in experiment or in disease, K and p are 
added together in equal quantities as the salt Kp. Whether p is a 
single species of anion does not matter as long as it is diffusible, uni- 
valent, and of uniform activity. The potassium K considered here is 
only that appearing with p beginning with the origin of the animal, 
or on mixing in the adult only that potassium associated with diffu- 
sible anions inside or outside the cell, that is, K = 'K,J + KE. Thus 
since p= K , p/W = K/W=K.. Substituting this into equation (2), 
with Cl and n held constant, we write 


eau, 
fa Tak 


Approximate values for the constants for the entire dog have 
been estimated from the literature. The lower curve of Figure 2 is 
drawn from Cl, = 31.4 K/(84.3 + K). 


(2) 


(3) 


Cl, 


(C1- M) (L+K)=(B-M)L 


FIGURE 2 
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. Further analysis reveals that the curve is an equilateral hyper- 
bola. If K > «, Cl, > the limit Cl, =CI=M.If K= —L, ChL=— o. 
From the derivative, 

d(Cl,) Resswod I, 


Cheeta) 2 ee 
the slopes at these respective values are 0 and w. The asymptotes, 
Cl, = M, K = —L are predicted (Figure 2). Visual inspection of the 
lower curve, on consideration of the conventional form L1Y1 = LoY2 = 
a?/2, gives 


(M-—Cl,)(L+K)=ML, (5) 


in which the left member results from substitution from the variable 
point P,, taking account of the fact that the asymptotes cross at 
(—L, M) rather than at the origin, and in which the right member 
results from substitution of the special point, the origin, which gives 
the constants. It is seen that a, (from the conventional equation) = 
V2ML. Conversion of equation (5) into equation (3) confirms the 
analysis. 

The Relation of Cl to K. Write Cl = (Cl/d)d. Then by substi- 
tution of d= N-+K, 


Cl= (Cl/d) (N + K). (6) 


In this form equation (2) is a meaningless expression of the equality 
between positive and negative radicals in extracellular fluid that holds 
in any theory. However, in the development to follow a useful rela- 
tion will appear between Cland Kk . 

Continuing the development as with equation (2), write 

Cl/W M(N + K) 
|= ——__ (N + K) = ——_—. (7) 
n/W + p/W Gers 

For this to follow Donnan relations, because of the substitution 
p/W =K, d, must again equal K as for equation (1). 

Equation (7) is an equilateral hyperbola with the same asymp- 
totes as equation (5), and from the upper curve of Figure 2, 


(Cl— M)(K+L)=(B-M)L; (8) 


a=\/2(B—M)L. Equation (8) converts to (7) on use of the identity 
BL = MA. While the two hyperbolas (5). and (8) have the same 
asymptotes, their foci are at different distances from the center. They 
are typical conjugates only when 2M = B, when FE =I, a condition 


not normally found in mammals. 
The Relation of Cell Chloride to Extracellular Chloride. Grant- 
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ing homogeneity of chloride in F and I, regardless of the mechanism 
of chloride transfer, one writes Cl EF + Cl,f=Cl. Whence, 

Ch = CUT E/LCl= D> L7i Cr: (9) 
or 

Cl=CUV/E —1/ECl=B—TI/E Cl. (10) 
These convert to either of the slope equations derived geometrically in 
Figure 3. Also dividing relation (5) by relation (8), we write 


(M—¢€L)/(Cl=M) = Mya): (11) 


FIGURE 3 


The Relation between the Three Variables. Division of equation 
(3) by equation (7) yields 


Cl,/Cl=K/(N+ K). (12) 


Variation of K with p Held Constant. If K varies without chang- 
ing p, the constancy of p/W eliminates K from the denominators of 
equations (3) and (7). Linearity results except for complicating 
volume changes avoided when K varies through addition of Kp. If 
K is substituted in part for N (Wilde, 1943), the hypotonic N allows 
I to swell at the expense of H. Suitable calculations show that (N + K) 
rather than N of expression (7) is a constant. The concentrations 
Cl,d,andp remain constant, while Cl, is a true linear function of K . 

Applications. Equations (3) and (7) may apply to entire ani- 
mals or to tissues excised after equilibration. Any anion supposedly 
diffusible may be substituted for Cl, for instance HCO;, SCN, Br, 
etc. Application to tissues, as muscle, is simpler since the constants 
L and M can be measured directly. Values of Cl, calculated from ex- 
pression (3) may be compared with Cl, estimated by some analogue of 
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inulin-space. Cl, calculated from equation (12) may also be so com- 
pared. In another approach, since the reciprocal of expression (3) 
is linear, 

(1/Cl,) =1/M + L/M(1/K), (18) 


the constants may be estimated by least squares. 

Estimation of the constants for equation (7) is more involved, 
especially if applied to an entire animal, for the behavior of chloride 
in different tissues may vary. However, the equation is very useful 
for predicting changes in plasma chloride in terms of plasma potas- 
sium when the latter increases in disease. For instance, from the 
equation plotted as the upper curve of Figure 2, Cl = 31.4 (120 + K)/ 
(34.38 + K), it can be calculated that if in such conditions as shock, 
hemorrhage, or renal disease the plasma potassium rises from 5 up 
to 10 milliequivalents per liter of water, presumably as Kp appears 
in the system, the plasma chloride may fall from 100 to 92 milliequi- 
valents per liter. No other explanation has been offered for such 
changes, which must represent a redistribution of chloride since none 
of the ion is lost to the exterior. 

If in the live animal only muscle participates in these relations, 
then I is reduced one-half and —ACI correspondingly. Also any ex- 
change between » and corpuscular Ci will reduce the change. On the 
other hand, dehydration of the E’ phase of the animal before injec- 
tion of the salt Kp , in that the L/I ratio is lowered, will increase the 
experimental drop in plasma chloride. 
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ON THE FORM AND STRENGTH OF TREES: 


PART I. THE TRUNK 
J. OPATOWSKI 
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_ The theory of strength of materials and known experimental rela- 
tions between the strength and the specific gravity of green wood are 
Vat to derive relations between some geometric characteristics of the 

ee. 


1. Some fundamental facts. N. Rashevsky (1943) has shown 
that a theory of form of plants can be based on concepts of metabo- 
lism and mechanical strength. It is the purpose of the present paper 
to indicate the contributions to be expected from the theory of the 
strength of materials. To derive numerical formulae, use is made of 
the results obtained by the Forest Products Laboratory of the U. S. 
Department of Agriculture on the elastic and strength properties of 
wood. These experiments are so exhaustive and cover such a wide 
variety of tree species that their use for theoretical purposes appears 
justified. The relations that we will need are the following ‘(Mark- 
wardt-Wilson, 1935, pp. 36-40): 

For compression parallel to the grain, the maximum strength o, 
is given by 


o- — 473G kg cm”, (1) 
and the modulus of elasticity #. is given by 
E,.= 204600G kg cm”. (2) 
For static bending, the maximum strength o is given by 
on = 1237G°“ kg cm”, (3) 
and the modulus of elasticity EH, is given by 
E, = 165900G kg cm”, (4) 


| where G is the specific gravity of oven dry wood based on the volume 
in green state. The units are kilograms and centimeters throughout 
the paper, unless otherwise stated. Equations (1) to (4) are results 
of tests obtained on green wood; they are applicable to live trees, at 
least for most species grown in this country, because these species do 
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not show differences in mechanical properties of sapwood, which con- 
sists mostly of living cells, and heartwood, containing dead cells only 
(Markwardt-Wilson, 1935, pp. 41,42). 

In calculating the weight of the tree we have to take into ac- 
count its humidity content. Therefore, we will put for the average 
specific weight of the tree, w = mG/1000 kg cm-*, where the factor m 
may vary between about 1.3 and 2.5 (Markwardt-Wilson, 1935, p. 7). 
For the actual compressive stress in the tree, when no external forces 
are acting, we shall write o./f, where f, a safety factor, is greater 
than one. 

2. A simple geometric model of a tree. Consider the trunk as a 
circular cone, and take into account the weight of the branches by 
multiplying the specific weight of the trunk by the ratio p = B/T, 
where B is the weight of all the branches and T the weight of the 
trunk (Greenhill, 1881, pp. 69, 70; Olsson, 1942, pp. 173, 174). If 
we treat the trunk as a simply compressed column (Timoshenko, 1940, 
p. 4), the calculation of the stress at the base of the trunk gives 
m(1 + p)fHG = 30 o, where H is the height of the tree in meters. 
Taking into account expression (1), we obtain from here 


m (1+ p) fH =14190. (5) 


Since f > > 1, equation (5) gives a limitation for the height of the 

tree which is seen to be independent of the specific gravity that is of 

the species of the tree as well as of the size of the base. To see the 

order of magnitude of f at its possible lower limit, consider a case 

of a very tall tree, a sequoia for instance: H = 100 meters, m = 2, 
=1. We find f = 35.5. 

For a column of the size of a tree the conditions of elastic sta- 
bility are likely to be more important than those of compressive 
strength (Timoshenko, 1936, pp. 64, 65; 1941, pp. 184-186). For a 
circular cone under the action of its own weight, the maximum height 
H, compatible with elastic stability is (Greenhill, 1881 pp. 69, 70): 
H, = (7.63E, r2/w)?, where r is the radius of the base. pee re- 
lation (4), we thus obtain 


[m(1 + p)/r?]?2 H,= 1082. 


One may call stability of the tree the ratio € = (H, — H)/H, where 
H is the actual height of the tree. For a very high tree, for instance a 
Sequoia sempervirens at its upper limit of size (Rehder, 1940, p. 48), 
H = 110 meters, r = 4 meters, so that by assuming m= 2, p=1, 
we get H, = 370, €= 0.7. 

3. The effect of eccentricity. One may take into account the ec- 
centricity of the tree with respect to its base by calculating the maxi- 
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mum compressive stress at the base due to the combined effect of com- 
pression and bending. We assume the base to be circular of radius 
r , and er the horizontal distance of the center of gravity of the tree 
from the center of the base. Then the maximum compressive stress 
at the base is (Timoshenko, 1940, pp. 230, 231): o=o6,(1 + 4 «), 
where o, is the compressive stress for a perfectly symmetric tree of 
the same weight and base. It is seen that the eccentricity increases 
the stress quite substantially. This is the reason why high trees are 
in general more symmetric. 

4. A model of a tree of maximum strength. Instead of assuming 
the tree is a cone, there is a possibility of deriving its shape from a 
general principle, according to which the material from which a plant 
is made is distributed in such a manner as to give a maximum 
strength. This principle was widely discussed by S. Schwendener 
(Holtermann, 1909; Fischer, 1904, pp. 119, 120) who gave actual evi- 
dence for its occurrence in the botanical world. We limit ourselves 
here to applying this principle to the trunk, which we treat as a com- 
pressed column. The principle of maximum strength is alone not suf- 
ficient to determine the shape. We add the condition of maximum 
height which is biologically justified at least in the case of closely 
grown trees. Since we keep out of consideration the phenomena of 
growth, we may assume the amount of material involved to be an 
assigned quantity. Then it is clear that the condition of maximum 
strength and maximum height implies that the stress be the same 
in each cross-section of the trunk, that is, that the trunk be of equal 
strength (Timoshenko, 1940, pp. 17, 18). 

We may take into account the weight of the branches in the same 
way as previously. Then if a(z) is the area of the cross-section of the 
trunk at a distance z from the base, the condition of equal strength 
gives: 

fq + p)w f¥a(z)dz= oa A(z). (6) 


The differentiation of this with respect to z gives dloga/dz = const. 
which may be easily integrated, but such an integral satisfies equation 
(6) only if the height H of the tree is infinity. This inacceptable con- 
clusion is a consequence of the assumption that-we made, that is, that 
the cross-section a is a differentiable function of the distance z. We 
have to abandon such conception of a smooth profile of the tree trunk 
and bring ourselves closer to the actual shape of the tree. 

Number the branches of the tree progressively starting, for in- 
stance, from the top: 1, 2,---,7,-+"- 7”. Assume that at the trunk 
base no branches exist. Call the 7-th trunk element that one which lies 
between the branches i and i + 1. Let each i-th trunk element be 
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bounded by two horizontal areas: a; at the top and a,’ at the bottom: 
Two consecutive trunk elements are considered in direct connection 
with each other, so that the areas a;’ and a@;,, are superimposed with 
a; internal to ai;,. Therefore, at the connection of two trunk ele- 
ments the cross-section of the trunk has a discontinuity. We imagine 
that at this discontinuity a branch is attached to the trunk. Let H; 
be the height of the i-th trunk element, 7; its weight, and B; the 
weight of the i-th branch. Put @;/@in = a;, t; = fTi/oc, 6: = fBi/oc, 
_bi/ti = pi, 1 + pit? = Ri, o:/pia = Ai, and with H; in meters, 


h, = fwH /o, = mf; /4730; (7) 


t; and b; have the dimensions of areas, h; is dimensionless. The ratio 
pi may be called the ramification, and a; the cross-sectional reduction 
of the 7-th element. We apply now the principle of equal strength to 
each trunk element. From the theory of strength of materials we have 
the following relations (Timoshenko, 1940, pp. 17, 18): 


a,—a;,,+6;, @ =a,exph,, t,=a,—-—a,, (8) 


the quantities with negative subscripts are here by definition = 0. It 
may be noted that our present tree model must have a branch on its 
top because otherwise we would get h, = o. The elimination of a,’ 
and a’; between equations (8) gives 


ti, — GQ; exph, eg ates Db; =a, (a exphi_, ; 


which divided side by side leads to a dimensionless relation between 
the heights of two adjoining trunk elements; 


exph; = R; ai) Aja exph;4 : (9) 
From relation (9) we have 
exph, — exph, = KR. — R, — Ai exph,. (10) 


If we assume that the height of the trunk elements increases from the 
top to the bottom, we have h, > A, , and from expression (10) R, > R, 
Or p2 < p,, that is, the ramification is smaller for the second than for 
the first trunk element. This agrees with the fact that the second 
trunk element must provide for a metabolic flow not only to its own 
branch 2, but also to the top branch 1. By this consideration of me- 
tabolic flow, it is justified to assume that all p; decrease from the top 
to the bottom. Since each side of expression (10) is > 0, we have 
also R, — R, > A, or a, < 1 — (p2/p,), that is, a, < 1 which is intui- 
tive too. 

As important as the mechanical strength is for the life of a plant, 
it is not the only factor affecting its shape. Therefore, equations (8) 
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must be completed by some mathematical expressions of metabolism 
and growth, which would make it possible to relate to each other the 
various trunk elements as well as the branches to the trunk. It is pos- 
sible, however, to derive from equation (9) a relation between certain 
geometric characteristics of the tree as a whole without any additional 
hypothesis. In general, h,’s are sufficiently small [cf. equation (7)] 
so that exph; = 1 + h;. This means the exponential profile of each 
trunk element is approximated by a suitable truncated cone. Using 
this approximation, and adding all the equations (9) to each other, 
we get 

ipa hh, 2 —1) 0° =n (11) 


where h = >,” h; represents the height of the whole tree; 
n/p* = Sa" 1/pi, AS =D Ai(1 + hi) /B"* (1 + hi), a = AP pt, 


and p* is harmonic average of all ramifications and may be called the 
ramification of the tree. A* is a weighted average of the coefficients 
A;,’s, therefore a* may be called the average cross-sectional reduction 
of the tree trunk. Since one may expect a* < 1, we get from relation 
(11) a condition for p* which is independent of n , that is (1 + p*) h> 
1 +h, , or by equation (7) ; 


p* > 4730(mHf)-1 + H, H — 1, (12) 


where H and H, are in meters. 

For instance, if H = 100, H, = 2, m = 2, f = 35, condition 
(12) is satisfied for any p*. If, however, H = 25 m, p* must be at least 
about 1.8, unless f > 35. Qualitatively the condition (12) is intui- 
tive: small trees (i.e. small H) which do not have heavy branches 
(small p*) have a high safety factor f. In expression (11), 1 is very 
large with respect to h, so that a* is close to 1, which agrees also 
with qualitative observations of the trees. 
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In continuation of previous studies of the mathematical biophysics 
of visual perception in relation to the aesthetic evaluation of visual 
patterns, an expression for the total intensity of excitation in a. dis- 
criminating center as a function of the intensity of the peripheral stim- 
ulus is derived. This expression is applied to the case of aesthetic judg- 
ments of similar polygons of different sizes. The theoretical conclusions 
are tested experimentally by use of standard psychological scaling meth- 
ods. The theoretical predictions are found to be in agreement with the 
experimental results. 


In a previous paper (Rashevsky, 1942) we pointed out the im- 
portance of the function J(S), introduced by A.S. Householder (1939) 
for the mathematical biophysics of perception. Previous notations 
(Rashevsky, 1940, pp. 153-154) are used in this paper. The quantity 
I(S), which is the total inhibition produced by cross-inhibitory fibers 
at any synapse, is proportional to the total net excitation 


E(S) = fo(s,h) dh, 


1 
at all the synapses, so that E(s) = 5 I(S), where 4 has the same 


meaning as before (Rashevsky, 1940). The function #(S) determines 
the total intensity E of the whole discriminatory center as a function 
of the intensity S of the peripheral stimulus. Using again the same 
notations as Householder (1939), we have the following expression 
for the function E(S): 


CS, h)dh ; 
E(S) = ToT (1) 
If, with Householder (1939), we consider the simple case 
(S,h) =h(S— A), (2) 
and define a quantity x by (Householder, 1939): 
St=S— 2h, = 2h, —'S, (3) 
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we have (Householder, 1939; Rashevsky, 1940) : 


Se or ee) 
he S, h dh = ’ (4) 
SaS 12 
and 
Hence 
Soe) 
LOSS (6) 
12(1+4S 2) 
The quantity « satisfies the equation: 
2A 
us? + «x? —1—0; pacer (7) 


Equations (6) and (7) give 


LE Suet) 9 
Bq eae re 
A 32 (1 +500) 


For large values of u, when u >> 1, I(u) varies as u?, since 
for such values of u, x is practically constant, decreasing slowly with 
u. This can be seen if we plot x as a function of wu, obtained by solv- 
ing equation (7) graphically, for different values of wu. Hence, defin- 


ing a as the average value of (3 — x?) /12 over the actually occurring 
range, 


3— a 
a= : (9) 
124 
and we have, introducing S instead of wu into expression (8) 
E(S) =as?. “ ={10) 


An actual graph shows that a relation of practically the same 
form as equation (10) holds also for small values of uw, so that equa- 
tion (10) may be used as a good approximation over a rather wide 
range of values of S. 

Relation (10) may be applied to different problems in percep- 
tion. In this paper we shall apply it to the problem of aesthetic judg- 
ments of geometrically similar patterns of different sizes. 

Consider a series of similar polygons which are similarly oriented 
with respect to the observer, and which each consist of » equal sides 
(Rashevsky, 1938, 1940, 1942). As far as the excitation of L-centers 
is concerned, contemplation of each polygon results in a p-fold exci- 
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tation of one L-center, the total excitation being p times the quantity 
E(S). As for the intensity S of the corresponding peripheral stim- 
ulus, we assume it to be proportional to the length / of the side of the 
polygon (Householder, 1940) for not too large visual angles. 
The total excitation of the L center will thus be proportional to 
. Since in a series of similar and similarly oriented polygons, the 
xcitation of the V, A , A , and possible symmetry centers remains the 
same, the total central excitation for each similar polygon should 
vary as 
Eig A +B (11) 


It must, however, be remarked that the above relation cannot 
hold for extremely small values of 1. For 1 = 0 the polygon shrinks 
to a point, and the excitations of the V, H, and A-centers vanish. 
What happens in this case neurobiophysically shall be discussed in 
another paper. For the present it is sufficient that relation (11) holds 
as long as we still perceive a polygon as such. 

Besides the excitation of the V, H, A, and possible symmetry 
centers, the term A in relation (11) also contains the excitation due 
to the perception of the background on which the polygon is drawn. 
A presentation of a white card without any drawing on it results 
itself in a finite cortical excitation. 

Since, according to previous publications, H is proportional to 
the intensity of pleasure (positive or negative) associated with the 
contemplation of a given pattern, equation (11) should represent the 
pleasantness values in a series of similar and similarly oriented poly- 
gons. By using some standard psychological method of scaling, it 
should be possible to verify relation (11) experimentally. Before pro- 
ceeding to this, we should, however, consider a few other aspects of 
the situation. 

All our considerations on the visual perception are limited to 
patterns of such a size which permits the detailed perception of the 
whole pattern without necessitating any movements of the head or 
even of the body as a whole. If a pattern is so large that it can be 
contemplated successfully only in parts by turning the head in differ- 
ent directions, the movement of the head sets up an additional pat- 
tern of proprioceptive stimulus, and those must be added to the regu- 
lar visual or oculomotor stimuli in evaluating the pleasantness or un- 
pleasantness of the contemplated pattern. 

By restricting ourselves to sufficiently small patterns, we shall 
have a limiting factor, namely the size of the card, or more generally 
of the uniform background upon which the pattern is presented. When 
the over-all size of the pattern is equal or almost equal to the size of 
the background so that some elements of the pattern are very close 
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to the “edge” of the background, the perception of that edge brings 
in additional elements which inhibit the centers involved in the per- 
ception of the pattern proper. In scanning the pattern, the eye may 
“jump” off into scanning the edge of the background. If the back- 
ground is a card, as is actually the case in most experiments, then 
the surrounding back of the card adds also to the inhibiting effect. 
Thus the aesthetic value of the pattern as such will decrease due to 
the effect of the “edge” of the background. The law of the decrease 
must be determined by an extension of the theory somewhat along 
the lines suggested in Section 5 of a previous paper (Rashevsky, 
1942). The theory in its present form does not yet predict the quan- 
titative aspects of this decrease. Qualitatively, however, the theory 
does require such a decrease. 

Hence we should expect the following situation to hold: for 
geometrically similar and similarly oriented polygons drawn on cards 
of the same size, the aesthetic value of the polygon should increase 
parabolically with the size of the polygons, when the polygons are suf- 
ficiently smaller than the cards. As the size of the polygons ap- 
proaches the size of the card, the aesthetic value should begin to 
decrease. 

If we have two series of similar and similarly oriented polygons 
identical in all respects except that one series is printed on larger 
cards than the other, then the deviation from the parabolic law in 
the series printed on smaller cards should occur at smaller sizes of 
the polygons. 

Finally, the following considerations must be made. Even a blank 
card represents a visual stimulus which produces some central exci- 
tation and exerts an inhibitory effect on the centers involved in the 
perception of the pattern proper. Let S, and S. be the corresponding 
stimulus intensities for the patterns and the card. The quantity S, 
in general may be expected to be an increasing function of the size 
of the card. Hence the value of EF, as given by equation (11) may be 
expected to slightly increase as the size of the card decreases. The 
amount of that increase will depend on the parameters of the cross 
inhibitory fibers. Thus with decreasing size of the card, we should 
expect an increase in the coefficent B in equation (11). As regards 
A, amore elaborate theory is needed to predict what would happen 
to it since, as remarked above, it contains itself the excitation due to 
the presentation of the blank card. 

To verify all the above conclusions, the three polygons, repre- 
sented on Figure 1, were used. Each polygon was printed in a series 
of six different sizes such that the length of each side were 3 mm; 
0-45 mm; 0-6 mm; 0-75 mm; 0-9 mm and 122 mm. Five such gse- 
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FIGURE 1 


ries were made for each polygon. In one series each polygon was 
printed on a card having the size 12:5 X 17-5 cm; the second series 
of polygons was printed on cards 10 X 14 cm large; the third, on 
cards 8-33 < 11-7 cm, the fourth, on cards 7-8 X 10-9 em, and the fifth, 
on cards 6-25 X 8.75 cm. Thus all cards are geometrically similar 
rectangles. Altogether there were 90 cards. 

A simpler shape of the card, namely a square, would have been 
preferable. Inasmuch, héwever, as this experiment was a by-product 
of another experiment made for a different purpose, in which differ- 
ent other patterns were made on rectangular oblong cards, the above 
specified shapes of cards were used. The cards with the polygons were 
presented to 50 subjects. Each subject was asked to arrange the cards 
into 11 groups in order of their aesthetic preference. Using this pro- 
cedure, it is possible to compute the psychophysical aesthetic ratings 
of each card with a polygon from the methed of paired comparisons. 
The computations, however, are extremely elaborate and time-con- 
suming when as many as 90 stimuli are used. Therefore, as a suffici- 
ently good approximation the average group number of a given de- 
sign was computed and used as a measure of the average aesthetic 
value. Figure 2 represents the results for polygons 1b. The other 
polygons gave similar results. 

First of all, we see that on the largest card the aesthetic rating 
increases parabolically with the size J of the side of the polygon almost 
up to the largest size. A decrease is noticeable only for the very larg- 
est polygon. As the size of the card decreases, the parabolic relation 
breaks down at smaller sizes of the polygons, as is predicted by the 
theory. Finally, there is a slight increase of the steepness of the para- 
bola with decreasing size of the card, although this effect is not too 
pronounced. It is more pronounced for polygons 1b. Thus the one 
quantitative and two qualitative predictions of the theory are veri- 
fied. The constant A in equation (11) increases slightly with decreas- 


ing size of the card. 
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A mathematical analysis of the absorption of an inert gas by a 
heterogeneous system of n phases, e.g., a limb consisting of n tissues, 
is presented. The total uptake of gas, ¢(t), up to time t is given in terms 
of arterial concentration, cardiac delivery, blocd volume, and the volume, 
permeability, and partition coefficient of each tissue. The theory pre- 
dicts. how the uptake curve should change in shape under a variety of 
physiological conditions, and how from the numerical values of the con- 
See ae values of certain tissue constants, e.g. permeabilities, may be 
obtained. 


Introduction. Many problems in the physiology of intact ani- 
mals involve the absorption by the tissues of some substance which 
has entered the blood elsewhere. While experimental work along 
such lines can be and has been done with great care, the resulting 
data are not of the type that can be readily interpreted. Consequent- 
ly,.a need has arisen for a theoretical treatment of this subject. The 
present note is intended to submit such a consideration. 

We will attempt here to treat the absorption of an inert gas in 
the blood by the tissues of a limb, and will use a vocabulary spe- 
cific to this case. However, the mathematical results obtained are in- 
dependent of vocabulary, and obviously could apply to anatomically 
different but functionally similar situations. 

For the purpose of the analysis, we regard the limb to be con- 
stituted as follows: There exists a blood chamber of volume V, cm* 
(the blood volume). Into this chamber there is an inflow of blood at 
the rate of R cm? sec; because blood and tissue fluid are practically 
incompressible, there is a liquid (blood + lymph) outflow also equal 
to R. The entering blood contains the gas in question at concentra- 
tion C gm cm-. The (spatially) average concentration of gas in the 
blood chamber is z» gm cm-*. In physiological contact with the blood 
chamber are 7 distinct tissues. The average concentration of gas in 


the 7-th tissue is 7; gm cm”. 


* The material in this article should be construed only as the personal opin- 
ion‘of the writers and not as representing the opinion of the Navy Department 
officially. 
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The matter of penetration constants of the tissues deserves some 
attention here, even though the general aspects are thoroughly treat- 
ed by Rashevsky (1938) to whose work the reader is referred. When 
the intracellular “solvent” is the same as the extracellular “solvent”, 
then we write by Fick’s Law , 


Penetration Rate x mee (a; Sy = hes eee (1) 
(gm cm-? sec’) mi) 


where K is a constant of proportionality, S; is the absorbing surface 
area in cm?, 6 is the thickness of the cell membrane, and p’; is the 
partition coefficient of the gas between water and “membrane sub- 
stance’. The constants, Kp’;/d are lumped together into a single h; 
em sec“, which is then the permeability of the membrane. On the 
other hand, when the “solvents” of the two sides are different, then 
for the penetration rate we obtain: 


K Sip; 


Penetration Rate = 


(14-220) an Si (Xo — ai vi), (2). 
where a; is the partition coefficient of the gas between the external 
medium and the intracellular solvent. For some of the tissues, e.g. 
adipose, which we will consider, it is necessary to use equation (2) 
rather than equation (1). It is perhaps superfluous to point out that 
in the first case, the steady state will be attained for x = x; , while 
in our case it is attained when x — aix; , and generally speaking, the 
steady state concentration within the tissue will be different from the 
steady state concentration in the blood. 

Derivation. The setting up of our differential equations follows 
easily from the principle of material balance, which holds for the 
blood chamber and each tissue independently. Thus: 


A(Vo Xo) ree’ 

———=RC-—R2x,-— dA, S;(% — a; 2) 
dt 1 

ad(vV; Xi) 


di =h; Si (%@ — a; 2%) $= 1,2, -0-n. 


The substitution, yo = C — a and y; = Ci/a; — x; (i=1, 2, --: n) 
renders the equations homogeneous and somewhat simplifies the solu- 
tion. We have then instead the set: 


dy i=n 
= Ryo + Shi Si (Yo — a: yi) (3) 
dt i=1 


as Vo 
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dy; 
Vi =F Si (Yo — a4 41) @=1,2,--72). (4) 


We proceed to solve the set by the method of undetermined coeffici- 
ents. The series: 
j=n 
YiH=> Ci; ert (5) 


j=0 


will be a solution of equations (3) and (4) provided we can find ap- 
propriate expressions for the C;; and the k;. This can be done as fol- 
lows. The solutions (5) must, by hypothesis, satisfy equations (3) 
and (4). identically, which means that after substitution is made the 
coefficients of e*' must vanish. Thus from putting relation (5) into 
relation (4) we obtain the solutions: 


Ci = Co; (2=1,2,---). (6) 


On the other hand, we have yet to impose boundary conditions on 
equation (5). The two cases which concern the physiologist are ab- 
sorption (when t=—0,y,—C, and y;—C/a;) and desaturation, (when 
t=0, y—0, and y; — 0). Since our personal interest is in absorp- 
tion, we will develop the equations for that case; however, it is clear 
that no mathematical difficulty is involved in the alternate conditions. 
For our case, then, when t = 0 relation (5) becomes: 


j=n 
DB Co =C 
j= 
and by (6), 
j=n Co; C 


s (aoe f= 1, 2,-0.0. (7) 
j=0 Gi h; SG a Vi k; Qi h; Pere 


Equation (7) is therefore a linear set of (n+1) equations which can 
be solved for the C.; by the usual methods. Then from equation (6) 
all the other C;; are obtainable. Determinantal solution for Co; gives: 
Co; peed O88 fos k, anor a) (V; V2 Sieys Vn) Ai; A>; oral Anj/ hj (4 On +++ On) 
(Iq thy «+ Tim) (Sy So-+* Sa) (io — Bey) Uy — By) — (8) 
(kn — k;) = Aj;C : 
where there is no factor, k; — k;, in the denominator, and where, 
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Si 
Ai; =a; h; Pape — k;, 
V; 


a function which we will have occasion to refer to later. Finally, put- 
ting equation (5) into equation (3), noting relations (6) and (8), 
we obtain the characteristic equation whose (n+1) roots are the k;: 


i=n Oj h,? S;? i=" 
V, k; yt (08 (9) 
i=1 Qj h; Som V; k; i=1 


By the use of standard approximation methods, all the roots of equa- 
tion (9) can be found once we have numerical values for Vo, R, a, 
h;, Vi, and S;. However, we may prove two useful facts about the 
roots of relation (9) without recourse to numerical methods. 

I. All the roots are positive. This follows because if any one of 
them, say z, were negative, we could write the left hand member of 
equation (9) as: 


HOV ee Viel +R+3S%S 
z)= —_—_____—_——_ | — PA oe (oa iSi 
in 0 A, S; + Vi2| ; | i=1 / 
The greatest value which the left-hand bracket could ever have would 
be, (for |z| = 0), © hi S;, whence the greatest value F(z) could have 
would be —(V.|z| + RF); thus F(z) could never be 0, and therefore, 
contrary to our hypothesis k; must always be positive. 
II. If then +1 roots are, in descending order, r™, r-)), -.. r© 


N) 
and if the products, a; hi ad are, in descending order, p™, p(™, 


-- p™, then these quantities fall into the sequence, r™ > p™ > ra 
> p@y...r) > p® > r©, For suppose that 7” and 7” are roots such 
that r” > r’. By hypothesis they both satisfy equation (9). On sub- 
tracting F'(7’) = 0 from F'(r”) = 0, we obtain, 


Gas +S Ry ey : 
0 moe avi Eee S, 3, ; 
(ats # ) (a i ) 
Vi V; 


It must therefore be that some of the terms in the summation are 
negative, i.e., there is at least one tissue for which, 


S 
fe OSE 
V 
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Now imagine the roots k; arranged in descending order, r™ > r(-») > 

- > 7. By the reasoning just given, there must be a product, 
ah S/V, between every two successive roots, but the number of such 
positions, n, is exactly equal to the number of tissues; thus if we 
adopt the notation, p™ > p’) >... > p® for the descending se- 
quence of products, we see that the two sequences must fall into the 
combined sequence, 7” > p™ >... > p™ > r, and our theorem is 
proved. 

From I and II it follows immediately that all the A,; are posi- 
tive, since all the k’s are positive and since, by II, there are, for any 
k; as many Aj;’s > 0 as there are (k; — k;) > 0, and as many 4ij,’s < 0 
as there are (k; — k;) <0. However, the coefficients, Ci;/C , which 
we may call A;;, need not be positive (see equation 6), nor must 
linear combinations of the A;; necessarily be positive. 

Next we turn to a consideration of the derivatives of k; with re- 
spect to V., R, and the constants of the s-th tissue, a,, h,, Vs and 
S,. These follow from equations (8) and (9). If by w; we denote the 
function: 


i 
>-03 


i=n Si 2 
Vet Sah ¥s( 7, Aix? 


1=1 t 


then we may write the various partial derivatives compactly as: 


ieee ees oy .0 (10) 
oR 
k; 
Se ie? da 20 (11) 
0(h; Ss) 
ok; yeaa 
= 0 bs ¥o( we! (12) 
0 Qs 8 
0 Vo 
Snake 
Oke; ==! Wj k; ORE ( ae: Vi Og < 0 . (14) 
av; Vs 


The signs of these derivatives are evident enough. In connection 
with their numerical magnitudes it may be noted that all the deriva- 
tives are proportional to @;, a quantity which diminishes with in- 
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creasing blood volume. The quantity, S/V, figures prominently in 
these functions, both by itself and where it appears in 4;;. Its pres- 
ence stresses the physiological importance of the shape in which the 
tissue is disposed with respect to its blood supply. 

The derivatives of the A;; are complicated expressions whose 
signs and magnitudes cannot be ascertained from our work so far; 
consequently we have omitted them from this section. 

We are now in a position to discuss the general solutions of the 
problem. Combining equations (5), (6), and (8) we obtain the con- 
centrations of the gas in the blood and tissues as: 


Lo = C {1 — (Ano e-* + Ag, CM + --- + Aon E-F*) } 


Uy = C 7 <2 (Ai e-hot He Aca echt | Ses = la eFnt) 


Oly 


(15) 


1 
Ly = C iz a (Ano eels An emt +... 4 Ax em)| 5 
Multiplying each concentration by the corresponding volume, we ob- 
tain the total amount of gas in the limb: 


Vi 4 


Vo ty + Vis ot Vue [Vo t 2 to te 
; Ol + On 


a (Aoo Vo ar Ax Vee Ano Vojyex 
—(AinVot Ag Vi Peet S, V,)e 
ae (Aon Vo 1 Ain Vi te ee = Axe V,)e" ! . 


To expedite further discussion, we will adopt certain symbols for the 
quantities in this equation: We shall denote by 


¢ = total amount of gas in the limb => Vix 


; : i=n V; 
6 = effective volume of the limb = V, + 3 — 
i= Oj 


Q; = the coefficient of the j-th exponential = ¥ A;; V; 


1=0 
j=n 
y = the variable sum = © Q; e*'*. 
j=0 
Our equation then becomes simply: 


o=C(O— Q,e%— Q et —... — Q, 6) = C(6—w). (16) 
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The solutions, ¢(¢), the total uptake, and x;(t) the uptake of the ¢-th 
tissue, are characterized by the following properties: 


(1) They contain a number of exponential terms, e-*s', equal to the 
number of tissues plus one. 


(2) These exponential terms are all decaying terms, ie., k; > 0. 
The curves, ¢(¢) and x;(t) both approach asymptotes, Cé and C/a; 
respectively. 

(3) The coefficients of e**, Q; and A;,;, depend in general upon all 
the quantities, V,, R, Vi, a;, h; S;, and also on the shapes in which 
the tissue masses are disposed. Nothing is asserted about the signs of 
these coefficients. 


(4) If the physical properties of the tissues are not radically differ- 
ent, then the mean value of the k’s should be roughly equal to the 
mean value of ahS/V for the tissues. (The exact relationship is 
Theorem II.) 


(5) The change in shape of 4(¢) or 4: (t) as Vo, R, Vi,a;,, or hi S; 
is changed is governed mainly by the changes induced in the k;’s; 
therefor by equations (10) to (14), the curve approaches its asymp- 
tote more rapidly with increases in Rk, h; S;, or ai , but it approaches 
it more slowly with increases in V, or V;. 

(6) Both ordinates $(t) and x;(t) are proportional to the delivery 
concentration, C. 

This concludes the purely mathematical considerations. The ap- 
plication to experimental data is considered elsewhere (Smith and 
Morales, 1944). 
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